Single-side-hydrogenated graphene: Density functional theory predictions 
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Hydrogenation has proven to be an effective tool to open the bandgap of graphene. In the present 
density functional study we demonstrate that single-side-hydrogenated graphene is a semiconductor 
with an indirect bandgap of 1.89 eV, in between the gapless graphene and wide bandgap graphane. 
We show that its electronic structure and lattice characteristics are substantially different from those 
of graphene, graphone, or graphane. The lattice parameter and C-C bond length are found to be 
lengthened by 15% of those of graphene. Our binding energy analysis confirms that such a single 
, sided hydrogenation leads to thermodynamically stable material. 
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Graphene is unarguably one of the most novel materials of modern time. Graphene has potential application in a 
wide range of fieldsji but the lack of bandgap is one obstacle for utilizing it in semiconductor based devices. There 
are a number of investigations on bandgap opening and tuning using various techniques. Giovannetti et a£ had 
predicted theoretically that the substrate might induce a small bandgap. However, experimentally substrate induced 
,_!^ , bandgap is an issue of controversy^— A more accepted approach is by creating finite size nanoribbons out of graphene. 

The electronic confinement in nanoribbons results in a substantial bandgap opening^— Controlling the width of a 
nanoribbon is a technologically challenging task, which puts a restriction on this method for precise tuning. A more 
• ■ promising method is hydrogenation of graphene. 

The interest in hydrogenation of graphene is twofold; it can be used to tune the bandgap and it can also be harnessed 
as a hydrogen storage material. In an important theoretical work, Duplock et a& first showed that upon adsorption 
I 1 of hydrogen on graphene a bandgap can be opened. Even a single hydrogen atom may lead to bandgap of 1.25 eV. 
Later, Boukhvalov et a£Q extensively studied the adsorption of hydrogen atoms on multiple layers of graphene. Their 
density functional study brought out the nature of the magnetism as well as stability of graphene in the presence of 
adsorbents. More recently, in an intriguing experiment Balog et a£i demonstrated that the hydrogenation can indeed 
be used to open the bandgap. Their results confirmed that patterned hydrogenation can be systematically applied 
to open up a bandgap. By a novel way of semi-hydrogenating graphene, Zhou et aS^ proposed a new ferromagnetic 
material, graphone, which exhibits a small indirect bandgap (0.46 eV). As we shall demonstrate in this paper, complete 
hydrogenation eventually leads to an indirect but a larger band gap. 

All of the above reports deal with partial hydrogenation of graphene. In a different approach, Sofo et in 2007 
jd. predicted the existence of completely hydrogenated graphene and named it graphane.^ It was two years later that 
\ Elias et experimentally confirmed the existence of the same, which sparked tremendous interest*^— Recently, 
Chandrachud et carried out a systematic investigation of hydrogenation of graphene. They have analyzed the 
transition from graphene to graphane as a function of hydrogenation and confirmed that the amount of hydrogenation 
i determines the bandgap. 

t-H ' Although hydrogenation of graphene has been under a series of investigations, the research addressed either partial 
hydrogenation (as in graphone) or complete but double sided hydrogenation (as in graphane). A recent report by 

. t-H ' Muniz and Maroudas^i considers the effects of hydrogen coverage on a single side of graphene plane. However, their 
study is based on the classical methods and does not address the electronic structure. At present, there is a compelling 
necessity to study in detail the electronic structure of single sided hydrogenation of graphene. It stems from the fact 

. . that when graphene is grown on a substrate like Si02 it is impossible for hydrogen to diffuse through the graphene 
to form double sided hydrogenatio n 14 ' 22 and the only possibility is partial or complete hydrogenation from a single 
side. Neglecting the substrate effects on graphene which were found to be relatively small^— we concentrate in this 
report on the electronic structure properties characterizing such single-side-hydrogenated graphene, which we shall 
call SSHGraphene. 

Our results bring out some remarkable properties of this material such as the indirect bandgap in its band structure. 
We also analyze the density of states (DOS) and their angular momentum decomposed counterparts along with the 
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orbital charge density, in order to exhibit the nature of its electronic structure. 

We carried out full unit cell optimization and detailed electronic structure calculations of SSHGraphene using 
density functional theory (DFT) with a plane wave basis as implemented in the Quantum ESPRESSO package^. We 
used ultrasoft pseudopotentials^ and the Generalized Gradient Approximation (GGA) with the exchange-correlation 
functional of Perdew, Burke, and Ernzerhof (PBE)i^ The plane wave basis cut-off was 60 Ry (~ 816 eV) for wave 
functions and 240 Ry (~ 3265 eV) for charge density, sufficiently large for graphene-like systems. A convergence 
threshold was set at 10 -8 eV for the energy self-consistency and 0.005 eV/A for the forces. Further, to maintain 
the accuracy, integration over the Brillouin zone was performed on regular 26 x 26 x 1 Monkhorst-Pack grids. The 
band structure was plotted on the lines joining the M, T, K, and M points, and the individual line segments were 
sampled using 50 grid points each. The corresponding precision was also maintained for the cell optimization carried 
out using the BFGS quasi-Newton algorithm. The convergence threshold on the pressure was kept at 0.1 kBar. The 
computational unit cell consisted of two carbons and two hydrogens. A vacuum space of 12 A was kept normal to 
the (SSHGraphene) plane to avoid any interactions between the adjacent sheets. We also validated our results for the 
bandgaps, bond lengths, and lattice parameters by making a comparison of different exchange-correlation functionals 
(see the Supplementary Material). 




FIG. 1: Top view of the SSHGraphene hexagonal structure (the solid line rhombus marks the unit cell) with carbon and 
hydrogen atoms shown in dark and lighter shade, respectively. A close-up with the bond lengths in Angstroms. 

We first discuss the geometry and then the electronic structure of single-side-hydrogenated graphene (SSH- 
Graphene). Whenever possible we shall put our results in perspective of graphene and graphane. For such a comparison 
it is worthwhile to review some properties of graphene and graphane before we discuss SSHGraphene. 

Graphene is a one-atom-thick sheet of sp 2 -bonded carbon atoms that are densely packed in a bipartite crystal lattice. 
It has two atoms per unit cell which has the lattice parameter of 2.46 A, with a carbon-carbon bond length of 1.42 
A. Although graphane is bipartite and hexagonal, its unit cell has four atoms (two carbons and two hydrogens) and 
has a larger lattice parameter, namely 2.51 AJ^ In graphane every alternate carbon atom is attached to a hydrogen 
atom from alternate sides of the plane. In response to the addition of hydrogens, the carbon atoms are displaced out 
of the plane towards hydrogen atoms. In short, the carbon atoms in graphane are no longer planar. 

The unit cell of SSHGraphene also contains four atoms, two carbons and two hydrogens. We carried out full 
optimization of the unit cell, including both the unit cell geometry and the atomic positions. The optimized geometry 
of SSHGraphene is shown in Figure Q] depicting a top view with the unit cell marked with the solid line and a close- 
up with the bond lengths. As seen from the figure the cell is similar to that of graphene, except that the lattice 
parameter for SSHGraphene is now enlarged to 2.83 A, which is larger than graphane (2.51 A) as well. Notice that 
the enhancement is necessary in order to accommodate the hydrogen atoms, as the unoptimized unit cell of graphene 
does not favor the complete hydrogenation. The increase in the lattice parameter is due to the increase in the carbon- 
carbon bonds. As seen from the zoomed-in view, the carbon-carbon bond length is increased from 1.42 A (in graphene) 
to 1.64 A. The lattice is thus expanded by 15%. The increase in the bond length upon hydrogenation is not surprising, 
as the same effect is also seen in graphane. Similarly, as expected, upon single sided hydrogenation the carbon atoms 
remain in one plane with the hydrogens forming another plane at 1.08 A. This is a typical bond length of C-H when 
bonded covalently. (In methane, for example, the bond length is 1.09 A.) To summarize, a comparison of (available) 
structural parameters of graphene, graphone, graphane, and SSHGraphene are given in table HI It also shows the 
binding energy per atom which is the signature of thermodynamic stability of the system. The biding energy for 
SSHGraphene is calculated using the pseudoatomic energies of carbon (Eq) and hydrogen (-Eh) atoms, and using 
AE = E c + En - -EsSHGraphene, where -EssHGraphcnc is the total energy of SSHGraphene. Thus, the higher the energy 
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a (A) 


2.46 




2.51 


2.83 


c-c (A) 


1.42 


1.50 


1.52 


1.64 


C-H (A) 




1.16 


1.11 


1.08 


AS/atom (eV) 


9.56 




6.56 


5.54 



TABLE I: A comparison of graphene and SSHGrapheneversus graphone and graphane as reported in the literature.— 1 ^ a is 
the lattice parameter, and AE is the binding energy. 
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FIG. 2: Band structure of SSHGraphene with the 1.89 eV indirect bandgap, versus graphane. Fermi level (dotted line) half 
way between the valance and conduction bands. The bold line marks the band which moved up in graphane, compared to 
SSHGraphene. 

the more stable the system. The binding energies for graphene and graphane are as reported in the literature.— The 
overall trend is quite straightforward. Graphene, having the smallest C-C bond, is the most stable of all. Although 
not as stable as others, SSHGraphene is still strongly bound. To put it in perspective, recall that benzene has the 
binding energy 6.49 eV/atom while acetylene has 5.90 eV/aton>H, and both are among the most stable hydrocarbons. 
Thus there is no doubt that SSHGraphene is indeed thermodynamically very stable. 

We now approach the electronic structure of SSHGraphene. It is well known that the graphene band structure is 
very sensitive to deformations of any kind. As noted before, there is a clear evidence that upon partial hydrogenation 
the bandgap of graphene is opened. It is thus easy to conjecture that the SSHGraphene would be a semiconductor. 
However the most remarkable feature of SSHGraphene is that it is a semiconductor with an indirect bandgap. At 
this point we recall that graphene is gapless, while graphane has a direct bandgap of 3.5 eV. 

The band structure of SSHGraphene shown in the upper part of Figure [2] clearly exhibits an indirect bandgap of 1.89 
eV. This is interesting, as this value lies in between the gapless graphene and the rather wide bandgap graphane (3.5 eV 
by DFT and 5.4 eV by GW method.—) Thus, SSHGraphene becomes a preferred organic candidate for semiconductor 
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FIG. 3: Angular momentum decomposed DOS of SSHGraphene and grapheme. p z from carbon and s from hydrogen contribute 
substantially more than p x and p y (identical to p x , hence not shown). Fermi energy is taken as zero. Gaussian broadening of 
0.1 eV is applied for better visualization. 



based devices. As mentioned earlier, a small indirect bandgap was also seen by Zhou et aS^ in graphone. However, 
we would like to point out that they did not report any optimization of the unit cell, which may considerably affect 
the bandgap as well as lattice parameters and bond lengths shown in Tabic [IJ 

The origin of the indirect bandgap can be inferred when we compare the band structures of SSHGraphene and 
graphane in Figure [5J A close examination reveals that the most significant difference consists in the position of 
one band of graphane marked in orange. In fact, this is the band that arises in graphene from the p z orbitals on 
carbon atoms and is degenerate at the K point. Upon hydrogenation this band lowers in energy and the nature of 
hydrogenation determines the lowering. 

In the case of graphane this band is significantly lower in energy as the hydrogenation is done in the best possible 
way, which is reflected in the stability seen in Table |TJ On the other hand, in the single sided hydrogenation all the 
hydrogen atoms are held at the same side, which naturally results in a repulsion among the hydrogen atoms. This 
further causes the crystal to expand. As a result of enhanced repulsion among the hydrogen atoms, the band marked 
in bold shifts upward in energy, which gives rise to an indirect bandgap. In short, the peculiar band structure is due 
to the enhanced repulsion among the hydrogen atoms which results in the shift of the most sensitive p z type band 
of the system. This also suggests that we may expect the p z type character associated with this band in the case 
of SSHGraphene as well. We verify the angular momentum character of the band by examining angular momentum 
decomposed density of states and by the partial charge density of the given band. Figure [3] presents the angular 
momentum decomposed DOS for SSHGraphene and graphane, showing the p x (identical to p y ) and p z components of 
DOS on the carbon atom, along with the s component on the hydrogen atom. It becomes evident from the comparison 
that in the case of SSHGraphene, the p z of carbon atoms and s of hydrogen atoms have shifted towards the Fermi 
level. In fact, the contribution from p z is almost twice that of p x . This observation strengthens our earlier conjecture 
that the band moved upwards in energy. 

The nature of DOS of SSHGraphene also underlines a strong hybridization of p z of carbon atoms with s of hydrogen 
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FIG. 4: Charge density isosurfaces of the band just below the Fermi level. The isosurface (at 0.1 of the maximum) covers the 
neighboring hydrogen atoms. Observe that the hydrogen atoms are connected along the arm-chair direction. 

atoms. This is the case, indeed, as seen from scrutinizing the individual bands. Figure 2] shows the charge density 
isosurface of the top most occupied band. The density is plotted at the T point. There are two clear observations to 
be drawn. First, there is a a bond established between the p z of carbon and s of hydrogen. Second, there is an s-s 
overlap between the two neighboring hydrogen atoms, established only along the arm-chair direction but not along 
the zig-zag. 

At this point, it is worth noting that it may be inappropriate to call this covalently bonded system sp 3 hybridized, 
since the characteristic bond angle of 109.5° is not present anywhere. Generally in the case of a few hydrogen atoms 
interacting with graphene or even for graphane, the underlying carbon atoms are displaced from their locations, 
showing the signature of sp 3 bonded system. However, in SSHGraphene all the carbon atoms remain in one plane, 
making it difficult to call it sp 3 hybridized. 

In summary, we investigated the single sided hydrogenation of graphene using first-principle DFT, and demon- 
strated that upon complete relaxation of the unit cell and atomic degrees of freedom SSHGraphene exhibits an 
indirect bandgap of 1.89 eV. This is appealing, as graphene is gapless while grapheme is rather wide bandgap 
semiconductor, making SSHGraphene a more desirable semiconductor. Our comparison of the band structures of 
SSHGraphene and graphane shows that the p z band located in graphane significantly below the Fermi level moves 
up in SSHGraphene, becoming the top most occupied band. This is attributed to an s-s overlap between the 
hydrogen atoms. Our analysis of the binding energies clearly indicates stability of SSHGraphene. Further accuracy 
check of the bandgap, bond lengths, and stability of SSHGraphene we have predicted here can be addressed when 
experimental data will be available. 

This work was supported by the Natural Sciences and Engineering Research Council (NSERC) and National Re- 
search Council (NRC) of Canada. The research has been enabled by the use of computing resources provided by 
WestGrid and Compute/Calcul Canada. B.S.P. also thanks D. G. Kanhere, K. Joshi and D. Saraf for their valuable 
inputs. Some figures are generated using VESTA. 27 
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A. Computational details 

All the calculations are based on DFT using a plane 
wave basis as implemented in the Quantum ESPRESSO 
package^. We use ultrasoft pseudopotentials 2 - to rep- 
resent the interaction between ionic cores and valence 
electrons. The Generalized Gradient Approximation 
(GGA) is used with the exchange-correlation functional 
of Perdew, Burke and Ernzerhof (PBE). 3 Kinetic energy 
cut-off on the plane wave basis for wave functions is set 
to be 60 Ry (~ 816 eV) and for charge density is 240 Ry 
(~ 3265 eV) which are sufficiently large for graphene- 
like systems. A strict convergence threshold for the self 
consistency is set at 10 ~ 8 while that for the forces is 
0.005 eV/A. Further, to maintain the accuracy, integra- 
tion over the Brillouin zone is performed on regular 26 x 
26 x 1 Monkhorst-Pack grids. The band structure is plot- 
ted on the lines joining M, T, K and M points, and the 
individual line segments are sampled using 50 grid points 
each. The higher level of precision is also maintained for 
the cell optimization which is carried out using the BFGS 
quasi-Newton algorithm with Hellmann-Feynman forces 
and the stress tensors. The convergence threshold on the 
pressure is kept at 0.1 kBar. The computational unit cell 
consists of two carbon and two hydrogen atoms. A vac- 
uum space of 12 A is kept normal to the plane to avoid 
any interactions between the adjacent sheets. 

To verify that the features seen in the band structure 
are not artifacts of the treatment on exchange-correlation 
functional, we have repeated the cell-optimization and 
band structure calculations using a different treatment 
to exchange correlation functionals including three differ- 
ent GGAs and one Local Density Approximation (LDA). 
Details of the results are given in the following section. 
Here we stress that our calculations confirmed that all are 
within reasonable agreement, except the bandgaps. This 
is understandable as the bandgaps obtained using DFT 
are sensitive to the choice of the functional. Since the 



qualitative picture remains unchanged, and since PBE is 
already proven to be reliable with graphene, grapheme 
and grapheme, it is appropriate to report the results ob- 
tained using GGA-PBE. 

B. Comparison of exchange-correlation functionals 

Choice of exchange-correlation functional plays a cru- 
cial role in DFT calculations. To test the accuracy 
of our result we now compare the results obtained us- 
ing different exchange correlation functionals. Except 
than the functional, other computational parameters re- 
main unchanged. We use three GGA based function- 
als, namely Perdew-Burke-Ernzerhof 3 , Heyd-Scuseria- 
Ernzerho£i£. and B3LYP 6 as well as one LDA based func- 
tional, Perdew-Zunger— for the comparison. The results 
obtained are summarized in table |TJ 

Because within DFT formalism the electron density is 
physically meaningful and the wave functions and energy 
bands are solutions of an eigenvalue problem for nonin- 
teracting electrons where their exchange and correlation 
interactions are shifted into the functionals, the bandgap 
becomes sensitive to the choice of functional. This is in- 
deed seen in the Table (|TJ . Although the bandgap values 
show considerable disagreement, other quantities remain 
in good agreement with each other. 

For a more conclusive value of the bandgap, treatment 
by a higher ab initio method is warranted. 





PBE 3 


HSE 4 ' 5 


B3LYR£ 


LDA 7 


Band Gap (eV) 


1.89 


1.35 


1.58 


1.22 


Lattice Parameter (A) 


2.84 


2.82 


2.87 


2.80 


C-C bond (A) 


1.638 


1.630 


1.656 


1.614 


C-H bond (A) 


1.076 


1.085 


1.070 


1.075 



TABLE I: Comparison of exchange-correlation functionals. 
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